Ghost imaging with a single detector 
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We experimentally demonstrate pseudothermal ghost imaging and ghost diffraction using only a 
single single-pixel detector. We achieve this by replacing the high resolution detector of the reference 
beam with a computation of the propagating field, following a recent proposal by Shapiro [J. H. 
Shapiro, arXiv:0807.2614 (2008)]. Since only a single detector is used, this provides an experimental 
evidence that pseudothermal ghost imaging does not rely on non-local quantum correlations. In 
addition, we show the depth-resolving capability of this ghost imaging technique. 



Ghost imaging (GI) has emerged a decade ago as an 
imaging technique which exploits the quantum nature 
of light. This field has attracted great interest, and has 
been in the focus of many studies since. In GI an ob- 
ject is imaged event hough the light which illuminates it 
is collected by a single-pixel detector which has no spa- 
tial resolution (a bucket detector) . This is done by using 
two spatially correlated beams. One of the beams illu- 
minates the object, and the photons transmitted by the 
object are collected by the bucket detector. The other 
beam impinges on a multipixel detector (e.g. a CCD 
camera), without ever passing through the object (the 
reference beam). Nevertheless, by correlating the inten- 
sities measured by the bucket detector with the inten- 
sities of each pixel in the multipixel detector, an image 
of the object is reconstructed [TJ. In a similar fashion, 
the diffraction pattern of the object can also be obtained 
('ghost diffraction', GD). In the first demonstrations of 
GI and GD the two beams were formed from a stream 
of entangled photons [2j [3]. The reconstruction of the 
image was attributed to the non-local quantum corre- 
lations between the photon pairs. Challenging this in- 
terpretation, Bennink et al. demonstrated GI using two 
classically correlated beams j3] , and triggered an ongoing 
effort to clarify the role of entanglement in GI and GD 
El U] HI H [10] [11] [Lj] [13] UJ]. It was soon discov- 
ered that many of the features obtained with entangled 
photons, are reproduced with a classical pseudothermal 
light source [9]. However, the nature of the spa- 
tial correlations exhibited with a pseudothermal source, 
and whether they can be interpreted as classical inten- 
sity correlations [T| [TTJ [13] or are fundamentally non-local 
quantum correlations [101 115) . is still under debate. 

In this Letter we experimentally study computational 
ghost imaging, a novel ghost imaging technique recently 
proposed by Shapiro [16]. In this technique the multi- 
pixel detector is replaced with a 'virtual detector', by 
calculating the propagation of the field of the reference 
beam. The image is reconstructed by correlating the cal- 
culated field patterns with the measured intensities at the 
object arm. Our measurements show that pseudothermal 
GI and GD can be performed with only one beam and 
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Figure 1: (Color online) Experimental setups for ghost imag- 
ing: (a) The standard pseudothermal two-detectors setup, 
where a ghost image of the object is obtained by correlat- 
ing the pseudothermal field measured by a CCD with the 
intensity measured by a bucket detector, (b) The computa- 
tional single-detector setup used in this work. A pseudother- 
mal light beam is generated by applying controllable phase 
masks tp r (x,y) with a spatial light modulator (SLM). The ob- 
ject image is obtained by correlating the intensity measured 
by the bucket detector, with the calculated field at the object 
plane. 

one detector. As noted by Shapiro, this proves that GI 
cannot possibly depend on any non-local quantum corre- 
lations. We also demonstrate scanningless 3D sectioning 
capability of this technique. 

Computational GI can be considered as a variant of 
the standard two-detectors pseudothermal GI. In pseu- 
dothermal GI, a spatially incoherent beam is generated 
by passing a laser beam through a rotating diffuser (Fig. 
[lji). The beam is then split on a beam-splitter, generat- 
ing the two spatially correlated beams required for GI. 
The essence of computational GI is to replace the rotat- 
ing diffuser with a computer controlled spatial light mod- 
ulator (SLM), which serves as a controlled phase-mask 
for the spatial phase of the light field ip{x,y) (Fig. [TJd) . 
A spatially incoherent beam is generated by applying 
pseudo-random phase patterns ip r (x,y) on the SLM. 
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Since for each phase realization r the controlled phase 
pattern is known, one can evaluate the field right after 
the SLM, E r (x,y,z = 0) = £(™) e *<Mz,y) (where E^ is 
the incident field on the SLM). Knowing E r (x, y, z = 0), 
the field at any distance z from the SLM can be computed 
using the Fresnel-Huygens propagator: 

E r (x,y,z) = j d^E^x-^y-^eW^, (1) 

where A is the wavelength of the source. In order to re- 
construct the transmission function of an object T(x, y) 
placed at z = L, the computed intensity patterns at the 
object plane I r = \E r (x,y,z = L)\ 2 are cross-correlated 
with the intensities measured by the bucket detector 
placed behind the object B r = J dxdyl r (x,y, L)T(x,y): 



G(x,y) = -^(B r -(B))I r (x,y) (2) 

r=l 

= (BI(x,y))-(B)(I(x,y)). 

Where (•) = A J2 r ' denotes an ensemble average over N 
phase realizations. Intuitively, one can see that the im- 
age is obtained by summing the calculated intensities I r 
with the appropriate weights B r : The larger the overlap 
between the generated intensity pattern and the trans- 
mission object, the higher is the intensity measured by 
the bucket detector B r , and thus the calculated I r (x,y) 
is summed with a larger weight. It is important to note 
that in conventional GI the object is also reconstructed 
according to Eq. The only difference is that here 

I r (x,y,z = L) are computed, whereas in conventional GI 
I r are obtained by measuring the intensities at the refer- 
ence arm using a multipixel detector placed at a specific 
distance z = L. Thus, similar to what is well known for 
conventional GI, one can show that G(x, y) is given by 
the transmission function of the object convolved with 
the coherence function of the field at the plane of the 
object [HUH]. The resolution of the reconstructed object 
is therefore limited by the coherence area of the field at 
the object plane. 

To demonstrate computational GI experimentally, we 
have constructed the setup presented in Fig. [IJd. The 
setup is based on a two-dimensional phase-only liquid 
crystal on silicon spatial-light-modulator (LCOS-SLM, 
Holoeye HEO1080P), with 1080x1920 addressable 8/imx 
pixels. The computer controlled 2D-SLM is illumi- 
nated by a CW Helium-Neon laser, which produces a 
Gaussian beam with a waist of wq = 740/im on the 2D- 
SLM plane. A 2cm x 2cm object (transmission plate) is 
placed at distance L = 84cm from the SLM, and a lens 
collects the transmitted light onto the bucket detector. 
In each realization a mask of 300 x 300 random phases is 
sent to the SLM, where each phase is realized by 3x3 SLM 
pixels. The reconstructed image using 16000 realizations 
is shown in Fig. |2ji, displaying the accurate reconstruc- 
tion of the object transmission T(x,y) (Fig. [2J1, inset). 




Figure 2: (Color online) Computational ghost image re- 
construction of a 2cm x 2cm transmission mask placed at 
L = 84cm. (a) Reconstructed image at the object plane, 
obtained with 16000 realizations. The inset shows the trans- 
mission mask, (b) A calculated intensity pattern of a single 
phase realization. The resolution of the reconstruction in (a) 
is dictated by the speckle size, (c) Reconstructed out-of-focus 
image, at a different z-plane (L = 15cm), demonstrating the 
depth-resolving capabilities of the computational method, (d) 
Measured SNR of the reconstructed image as a function of the 
number of realizations (blue dots). The theoretical line de- 
picts \/~N dependence. 



This experimental result cannot be attributed to nonlo- 
cal quantum correlations, since only a single detector was 
used. We note that one can reconstruct an object placed 
at any distance L from the source, both in the near- and 
far-field zones, as long as the field at the object plane can 
be calculated using Eq. (JlJ. This is demonstrated below, 
with an object placed at L — 11cm. 

An intensity pattern calculated from a single realiza- 
tion is shown in Fig. [2L, revealing the speckle field that 
impinges on the object for this specific realization of ran- 
dom phases |T7]- Reconstruction of an image at a z-plane 
different than the actual location of the object, results in 
an out-of-focus image of the object, indicating the depth- 
resolving capabilities of computational GI (Fig. |2p). This 
is analogous to scanning the z-position of the multipixel 
detector in conventional GI, but is performed here for 
all z-locations simultaneously, without the need for addi- 
tional measurements. This feature allows for a full three- 
dimensional reconstruction of the object field, with the 
only price of executing more calculations. 

Since the computational GI scheme is completely anal- 
ogous to pseudothermal GI, the two techniques share 
similar resolution and signal-to-noise (SNR) properties 
|16l [T71 [T8l [19]. To summarize these, the transverse res- 
olution is determined by the coherence length (speckle 
size) at the object plane, which is given by the Van- 
Cittert Zernike theorem Sx(z) = Xz/nw . The depth res- 
olution is given by Sz(z) — 2n8x 2 /X. For the presented 
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Figure 3: (Color online) Simultaneous reconstruction of an 
object and its diffraction pattern, with 8000 realizations, (a) 
Reconstructed ghost image of a two-slit transmission plate. 
The resolution of the image is determined by the speckle 
size, 5x ~ 0.018mm (b) Transmission microscope image of 
the double-slit slide. The width of each slit is 170/im and 
the separation is 400/mi. (c) Ghost diffraction pattern re- 
constructed using a pinhole detector. The resolution of the 
diffraction image is Sk ~ 1.4mm" 1 , lower than the Fourier 
limit 8x8k < 0.5. (d) A cross-section of the ghost diffraction 
pattern shown in (c) (dashed blue), and the theoretical two- 
slit diffraction pattern calculated from the dimensions of the 
double-slit (red). 

experimental parameters, Sx(L = 84cm) = 230/im and 
8z(L = 84cm) = 50cm, in agreement with the exper- 
imental results. In pseudothermal GI, the SNR scales 
as the square-root of the ratio of the number of realiza- 
tions N and the average number of speckles transmitted 
through the object N s , SNR cx yjN/N s [18] [19]. Since 
N s is given by the ratio of the object transmissive area to 
the coherence area, there is a clear trade-off between res- 
olution and SNR. The SNR as a function of the number 
of realizations for the reconstructed image of Fig. [2^,, is 
presented in Fig. [2ji, in agreement with the theoretical 
prediction. A movie visualizing the image build-up can 
found at [20] . 

One of the features in common to ghost imaging with 
a nonclassical source and a pseudothermal source, is 
the ability to resolve both the object and its diffrac- 
tion pattern with high resolution. In conventional GD, 
the diffraction pattern of the object is reconstructed by 
replacing the bucket detector with a small single-pixel 
('pinhole') detector, placed at the Fourier plane of the 
collecting lens [6] |8]. The GD is reconstructed by cor- 
relating the intensities measured by the pinhole detector 
with the diffraction pattern of the field in the reference 
arm. In fact, both the near-field object image (GI) and 
its diffraction pattern (GD) can be obtained by chang- 
ing only the optical setup at the reference arm. In our 
single detector configuration, since the reference arm is 
virtual, the only required change is in the computational 
procedure. One can therefore perform GI and GD si- 
multaneously, with a single set of measured data (i.e. 



the intensities measured with the pinhole detector B r ). 
Both images are reconstructed using Eq. where for 
GI I r is computed according to Eq. (JTJ, and for GD 
is obtained by calculating the Fourier transform (FT) of 
E r (x,y,z = 0), I r = \FT(E r (x,y,z = 0))| 2 . The intensi- 
ties measured by a pinhole detector placed on the optical 
axis are given by B r = | J dxdyE r (x,y, L)T(x,y)\ 2 . We 
note that by using a pinhole detector instead of a bucket 
detector, the SNR of the GI is degraded, since only a 
small fraction of the transmitted light is collected. 

In order to demonstrate the simultaneous reconstruc- 
tion of GI and GD, we have placed a double-slit trans- 
parency at L = Warn from the SLM and a small pinhole 
detector at the Fourier plane of the collecting lens (a 
single 8.6fim x 8.4/xm pixel of a CCD camera). We re- 
constructed the double-slit diffraction pattern using the 
pinhole-detector measured intensities, and its transmis- 
sion image using a bucket-detector measured intensities 
(summing over the CCD's pixels), obtained simultane- 
ously with the same phase realizations. We have also 
used the pinhole detector to reconstruct both the GD 
and GI, but with much lower SNR for the latter (not 
shown). Fig. [3] summarizes the results of the near-field 
and far-field image reconstructions for the double slit 
transparency, obtained with 8000 realizations. Compar- 
ing Fig. [3ji and Fig. [3]}, it is obvious that the product of 
the resolution of the near-field image (8x) and of the far- 
field image (Sk ~ l/^o) is much smaller than the Fourier 
limit, i.e. 8x8k — 0.025 -C 0.5. This cannot be obtained 
by simply scanning the object with a coherent laser beam 
j6]. However, it does not violate the Einstein-Podolsky- 
Rosen (EPR) bound for classical light, as was previously 
shown with a pseudothermal source |S]. Knowing both 
the near- and far-field intensity patterns with high resolu- 
tion enables the reconstruction of the object's transmis- 
sion amplitudes and phases, using the Gerchberg Saxton 
phase- retrieval algorithm [21]. This ability of computa- 
tional GI might be appealing for phase sensitive applica- 
tions. 

In conclusion, we have shown that pseudothermal GI 
can be performed with only a single detector, thus prov- 
ing that it does not rely on non-local quantum correla- 
tions. In addition, we have demonstrated the depth re- 
solving and diffraction imaging capabilities of this single- 
detector technique. Finally, we note that in computa- 
tional GI the main complexity is shifted from the exper- 
imental apparatus to the calculation. Thus, it allows for 
3D reconstruction by only post processing the retrieved 
data, eliminating mechanical scans. One might consider 
applying computational GI for other 3D imaging tasks 
which are not necessarily transmission-based. Two ex- 
amples for such are radar applications, and scanningless 
depth- resolved microscopy using fluorescent probes. Fur- 
thermore, the use of a SLM enables the implementation 
of closed-loop feedback schemes, potentially reducing the 
number of realizations needed for image reconstruction. 
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